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ABSTRACT 

A literature compilation of nuclear star cluster (NSC) masses is used to study the correlation 
between global and nuclear properties. A comparison of observational data to the predictions of semi- 
analytical galaxy formation models places constraints on the co-evolution of NSCs, massive black 
holes (MBHs) and host galaxies. Both data and theoretical predictions show an increased scatter in 
the NSC scaling correlations at high galaxy masses, and we show that this is due to the progressively 
more efficient ejection of stars from NSCs caused by MBH binaries in more massive stellar spheroids. 

Our results provide a natural explanation of why in nucleated galaxies hosting a MBH, the ratio 
(AInsc + AfMBH)/Afbuige (with Mbuige the host spheroid’s mass) shows significantly less scatter than 
.^Nsc/Afbuige, and suggest that the formation of MBHs and NSCs are not mutually exclusive, as also 
supported by observations of co-existing systems. Both MBHs and NSCs represent generic products of 
galaxy formation, with NSCs being destroyed or modified by the merger evolution of their companion 
MBHs. 

Subject headings: galaxies: Milky Way Galaxy - Nuclear Glusters - stellar dynamics 


1. INTRODUCTION 

Observations over the last two decades have un¬ 
veiled the existence of compact stellar nuclei at 
the centers o f most intermediate and low luminos¬ 
ity galaxies (CTarollo et a l.l Il998 t iBoker et al.l 120021 : 
Iden Brok et al.l I20l4 iGeorgiev fc Bokeil 120141 1. With 
masses of order 0.1% the stellar mass of the host galaxy 
and typical radii of only a few parsecs, these clusters, 
often referred to as NSGs, are the densest stella.r sys¬ 
tems observed in the Universe (e.g., lWalcher et al.l[2005l : 
iCote et nil2006l : iCarson et al.ll201,^ . 

The masses of NSCs (Mnsc) obey fairly tight 
correlations with properties of the host galaxies 
such as the bulge velocity di spersion (a) and mass 
(Afbuige) (|Ferrarese et al.l[2006all . A number of authors 
have shown that these correlations are much sh a llower 
than the corresponding ones for MBHs. iGrahaml (j2012f ) 
finds that NSC masses obey Mnsc ~ while the 
masses of MBHs, Mmbh, follow a much steeper scal- 
ing relation Mmbh tr® (|Ferrarese fc FordI 120051) . It 
is intriguing that for galaxies containing both a NSC 
and a MBH, the ratio Mcmo /Mbuige (with Mcmo = 
Mmbh + Mnsc the total mass in central massive ob¬ 
jects) shows le ss scatter than either Mmbh /Mbuige or 
MNSc/Mbuige (|Kormendv fc Hdl2013D . The facts that 
both NSCs and MBHs are found to follow tight corre¬ 
lations with their host galaxy properties, that they are 
found to coexist in some galaxies, and that in these galax¬ 
ies they have comparable masses, point toward a picture 
where NSCs and MBHs are both generic by-products of 
galaxy formation, and the growth mechanisms of NSCs 
and MBHs are related to one another and to the host 
galaxy evolution. 


We have developed a semi-analytical model following 
the formation and evolution of galaxies, MBHs and 
NSCs along cosmic history. Previous calculations 
assumed NSC formation to take place in isolated galax¬ 
ies (e.g., lAgywal fc Milosavlievid[20lH lAntoninil 120131 


Gnedin et al.! ^14 Arca-Sedda fc Capuzzo-Dolcettal 

2014 lAharon fc Peretsl 1201511 . thus neglecting the 
possible role of galaxy mergers, as well as in-situ star 
formation processes. Also, these early idealized attempts 
did not explore the possible interplay between MBH and 
NSC evolution. Our model sheds light on exactly these 
points, i.e. it allows us to assess for the first time the 
role of galaxy mergers, MBH mergers and nuclear star 
formation on the formation and evolution of NSCs. 


2. METHOD AND MODEL 

The backbone of this st udy is the sem i-analytical 
galaxy formation model of iBarajass Jwith the 

improyements described in iSesana et al.l ( ^14) ). This 
model tracks the eyolution of baryonic structures (both 
in central galaxies and in satellites) along dark matter 
merger trees produced with an extended Press-Schechter 
formalism, modified to reproduc e the results of A^-bod y 
simulations of dark matter halos (jParkinson et al.ll2008fl . 
The baryonic structures ultimately originate from the 
hot, largely unprocessed intergalactic medium, which 
cools to form cold gaseous galactic disks. These un¬ 
dergo quiescent star formation, which chemically enriches 
the interstellar medium and gives rise to stellar disks. 
Both stellar and gaseous disks are disrupted by bar in¬ 
stabilities - giving rise to pseudo-bulges that undergo 
a quiescent disk-like star-formation process - and by 
major mergers - which form classical bulges with vio¬ 
lent star-formation bursts, possibly triggered by turbu- 
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Fig. 1. — Comparison of our model to the observed NSC scaling relations. The upper and lower panels present respectively the predicted 
^NSC 3-iid Mnsc + ^MBH for the “erosion” model, which accounts for the mass deficit caused by MBH binaries, while the middle panels 
give Mnsc for tho “preservation” model, where this effect is absent. Blue lines are the median correlations predicted by our model; 
dot-dashed lines represent regions containing 70% and 90% of the NSCs produced in our models at a given stellar mass. Arrows are 
upper limits to the mass of NSCs. Red symbols denote galaxies with observationally constrained MBH masses; the red symbols in the 
different panels correspond to the same suloset of galaxies, but in the upper panels they give MnsCi in the bottom and middle panels they 
give Mnsc + ^ mbh - For galaxies with only an upper limit to MBH mass, we set Mmbh = Oj but star symbols show how these points 
would move if we were to assume MBH masses equal to these upper limits. The red line s are fits to the total N SC-fMBH masses, while 
dot-dashed-purple lines show for comparison the MBH — host galaxy scaling correlations (IKormendv fc H^I201.?li . 


lent gas flows. Besides further enriching the interstellar 
medium, bar instabilities and (mainly) star-bursts are 
assumed to funnel cold, chemically enriched gas into a 
low-angular momentum nuclear re servoir, available for 
accretion onto the central MBH (jGranato et al.l 12004 
iHaiman et al.l 12004) . Possible mechanisms capable of 
removing angular momentum from the gas are shocks 
(caused e.g. by merger s ) or radiation drag f rom star- 
bursts (jUmemnral l2004 ICranato et al.l 12004) , and we 
thus assume a linear correlation between the feeding of 
the nuclear reservoir and the star-formation rate in the 
bulge (or pseudo-bulge). The MBHs may form from ei¬ 
ther high-redshift seeds arising from Pop III stars (i.e. 
^ IOOMq “light” seeds; iMadau fc ReesI 120011) or from 


the remnants of protogalactic disks ( i .e. ^ lO^Mg) 

“heav y” seeds: liToushiaDDas et al.ir2004 iVolonteri et al.1 
120081) . In the former case, we allow moderately super- 
Eddington accretion onto the massive MBHs so as to 
better re produce the quasar luminosity function at high 
redshift (IMadau et al.ll20l3) . 

Here, we improve on this model by supplementing it 
with semi-analytical prescriptions for the formation and 
evolution of NSCs. These can grow ei ther by local sta r 
formation in the nuclear regio ns (e.g..lSeth et al.||2( I06D . 
or by infall of stellar clusters (jTremaine et al.lll975[ ) to- 
ward the center of the galaxy, as a result of dynamical 
friction. 

In more detail, we assume that the nuclear reservoir 














































3 


The imprint of black-hole mergers on nuclear star clusters 


from which the MBH accretes also undergoes quiescent 
disk-like star formation - for which we adopt the same 
star-formation prescri ption as in galact i c disks and pseu¬ 
dobulges, i.e. that of iKrumholz et id] (j2009tl - and we 
then assume that the stars formed in this nuclear reser¬ 
voir contribute to the NSC. Note that even though we 
model star formation in this nuclear gas reservoir as qui¬ 
escent, our models also assume that the reservoir’s feed¬ 
ing is correlated with star formation in the bulge, and 
thus mainly with turbulence-driven starbursts following 
major mergers. 

As for the infall of stellar clusters, we assume that the 
latter form with an efficiency /gc during star-formation 
episodes in both disks and bulges/pseudo-b ulges, and 
with initial mass function dn/dM oc M~^ (iBik et al. I 
1200311 ■ between M = 10^ Mq and M = 10^ Mq. Ini¬ 
tially, these stellar clusters follow the same spatial distri¬ 
bution as the stars, but after their formation we follow 
their (average) infall towards the NSC under dynamical 
friction from both the stars and gas. In doing so, we 
also account for both tidal heating and evaporation of 
the stellar clusters due to the galaxy’s tidal field, and for 
the tidal disruption by the MBH when they approach the 
center. When the MBH mass becomes > 10®“® Mq, the 
stellar clusters are completely disrupted before reaching 
the center, and the contribution of stellar cluster infall 
to NSC growth becomes negligible (IAntoninill2013[l . 

Whenever two galaxies merge, their MBHs are as¬ 
sumed to form a binary. Due to the interactions with 
the stars in the nuclear region, a MBH binary ejects 
from t he system a m ass (in stars) comparable to its own 
mass (|Merrittll2006f) . Also, when a MBH binary finally 
merges, the anisotropic emission of gravitational waves 
imparts a kick (o f up t o several thousands of km/s, cf. 
ICampanelli et all (I2007f) ') to the MBH remnant, and this 
also removes stars from the nuclear region. We thus 
model the total mass deficit Mej caused by the inspiral 
and merger of MBH binaries onto the NSC mass by 

Mej « 0.7gO ®Mbi„ + O.SMbin In ( — 

yOgr 

+5Mbi„(Wick/Ksc)^-^'' , (1) 

where g < 1 is the binary’s mass ratio, Mbin its mass. 
Oh the semi-major axis when the binary first becomes 
“hard” (i.e. tightly bound), Ogr the separation at which 
gravitational-wave emission starts driving the binary’s 
evolution (cf. explicit expressions for Oh and Ogr in 
iMerritd (|2013ll '). t4ick the post-merger kick velocity (com¬ 
puted using the fit to numerical- relativity simulations 
presented in ivan Meter et all (I20I0I1 1. and I4sc the escape 
velocity from the central parts of the galaxy (computed 
using the bulge and NSC density profiles). 

The first term in equation m accounts for t he mass 
defici t before the MBH binary becomes hard (iMerrittI 
1200611 ; the second is the mass ejected from Oh to 
Ogri and the third represents the mass deficit gen- 
erate d as the kicked MBH hea ts up the surrounding 
core (IGualandris &: Merrittl[2008ll . 

We also account for the tidal disruption of NSCs by 
MBHs during galaxy mergers. More precisely, in a 
merger between a galaxy containing a MBH and one con¬ 
taining only a NSC, the NSC is tidally disrupted before 


CO 



redshift 


Fig. 2. — Main progenitor evolutionary track of a galaxy whose 
NSC is partially disrupted by MBH mergers. We show the mass of 
the NSC, of the MBH, of the central gas reservoir, and the total 
mass in stellar clusters, as a function of z. 

reaching the center of the newly formed galaxy, if the 
MBH mass is > 10®“® Mq. This is the same process 
described above for the infall of stellar clusters. 

Next, we present our model’s predictions, with the free 
para meter f^r. set to its Milky Way value, i.e. /gc ~ 
0.07 (iKruiissenI 120121 1. We have checked the robust¬ 
ness of our results against the MBH seed model (i.e. 
light vs heavy seeds, with several halo occupation num¬ 
bers), a different value for /gc < 0.2, a variable /gc (set 
to 0.07, 0.04 and 0. 5 in disk, quiesc ent and starburst 
galaxies respectively: IKruiissenI [20121) and other details 
of our model (i.e. merger-tree resolution, initial redshift 
of the simulations, AGN feedback prescriptions, etc.). 
Indeed, we will show that the crucial ingredient to re¬ 
produce the observed global to NSC correlations is the 
scouring effect due to MBH binaries described by equa¬ 
tion O- More details of our mo del will be presented 
elsewhere (lAntonini et ahllin prep.ll . 


3. RESULTS 

We describe the results from two realizations of our 
galaxy-formation model, namely our fiducial model, 
which accounts for the “NSC erosion” due to MBH bi¬ 
naries through equation O; and a “NSC preservation” 
model where that effect is absent (i.e. we set Mej = 0 
throughout the entire cosmic history). 

Our model’s results are compared to a literature com¬ 
pilation of NSC masses. We constructed our sample 
of NSC objects by com bining data from Table 1 of 
IScott fc GrahamI (120131). Table 2 of lErwin fc Gadottil 
(j2012ll . Table 2 of iNeumaver fc Walched ( 20121 1 and bv 
estimating NSC masses for the early typ e galax ies in 
the Fornax cluster cat alog of [Turner et al.l (|2012f) . The 
masses of the NSCs in [Turner et al.l ( 201211 sample were 


calculated from the ( q — z) colors given in that paper and 
from the relations of IBell et al.l (1200311 . The galaxy ve¬ 
locity dispersion and bulge mass were also obtained from 
these pap ers or, when not ava ilable, from the Hyperleda 
database (|Paturel et al.ll200ll . 

The upper panels of Figure [T] compare the observed 
Mnsc — O' and Mnsc — Afguige relations to our “NSC ero- 
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Fig. 3. — Redshift evolution of NSC and MBH M — cr relations in our models with (left panels) and without (right panels) mass ejection 
due to MBH binaries. The bottom right panels compare our predictions to data. Blue open circles are NSCs; black open circles are MBH 
masses from ITremaine et al.l II2002I1 . Note the different evolution of the NSC mass vs cr relation in the two models, at z < 2. Since reliable 
MBH mass measurements have been obtained almost exclusively for main early type galaxies fe.g.. I^rrarese fc Fordir2005ll . in this analysis 
we have neglected MBHs and NSCs in satellite galaxies, and only included galaxies with bulge to total mass ratio B/T > 0.7, in order to 
compare with the observed MBH M — a relation. However, if all galaxies were included, the NSC scaling correlations would be very similar 
to those shown here (see Fig.[T] which was produced by considering all galaxies). 


sion” model. We note the excellent agreement of our pre¬ 
dictions with observations. In particular, the predicted 
scaling relations appear to broaden at higher velocity dis¬ 
persion, as do the data. The broadening of the data ap¬ 
pears to be caused by the subsample of galaxies with 
determined MBH mass (red points in Figure [T]). The 
nuclear to global correlations in these galaxies appear to 
flatten or even decline for spheroids with higher veloc¬ 
ity dispersion, suggesting that the MBHs have played an 
important role in affecting the properties of their host 
NSCs in these galaxies. 

The blue lines in the middle panels of Figure [T] shows 
the scaling relations obtained from the “NSC preser¬ 
vation” model. Note that the broadening at high it’s 
present in the upper panels disappears. This demon¬ 
strates that the broadening of the scaling relations seen 
in the “NSC erosion” model is a consequence of dynam¬ 
ical ejection of stars from coalescing MBH binaries. To¬ 
gether with the model’s predictions, we plot Mcmo = 
ATnsc + Mush vs a and Mbuige, for nucleated galax¬ 
ies with measured MBH mass or an upper limit to it. 
MBH masses and upper limits wer e taken from Table 1 of 
lErwin fc Gadottil (I2012D. table 2 ofiNeumaver fc Walche^ 
(j2012D . and Table 1 of iGraham fc Soitlerl (j2009t) . The 
lines show fits to these data. As can be seen, when the 


disruptive effect of MBH binaries is not included, our 
model recovers the fitted functional form of the relation 
between Mcmo and galaxy properties, both in slope and 
normalization. 

The fact that our “NSC preservation” model predicts 
scaling correlations that are consistent with relations in¬ 
volving Mcmo is not coincidental, but rather quite re¬ 
vealing. Indeed, the mass ejected by a MBH binary is 
of the order of the mass of the binary itself, with only a 
weak dependence on the binary’s mass ratio and the in i- 
tial density distribution of stars (jMerritt fc Szelll 120061) . 
Hence, the mass ejected in N MBH mergers is approxi¬ 
mately Mej « A^Mmbh, with Mmbh the hnal MBH mass. 

Because infall of gas to the nuclear region can 
rapidly rebuild a NSC, the relevant quantity to esti¬ 
mate is the number of MBH mergers since the era at 
which most of the gas was d epleted from the galaxy. 
iHaehnelt fc Kanffmannl (j2002[ ) compute this quantity us¬ 
ing semi-analytic models for galaxy mergers similar to 
ours. They find that the number of comparable-mass 
MBH mergers is weakly dependent on galaxy luminosity, 
and has a small dispersion around a median of 1, with 
values 4 even for the most luminous galaxies. Accord¬ 
ingly, derived mass deficits in “core galaxies” are also 
found to be peaked around a value of Ri Mmbh ijCrahaml 
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120041 . The total amount of ejected stellar mass is there¬ 
fore determined mainly by, and is roughly equal to, 
A^mbh- Therefore, Mcmo = -M^nsc + -^mbh can be re¬ 
garded as an approximate estimate (with uncertainties 
of order a few) of the NSC mass that was in place before 
mass erosion due to MBH binaries became important. 
Our models do indeed confirm this. 

These facts explain why the “NSC preservation” model 
predicts scaling correlations that are consistent with re¬ 
lations involving McmO: and provide further support for 
the idea that the broadening/bending of the observed 
correlations between NSC and host-galaxy properties is a 
consequence of the nuclear mass erosion caused by MBH 
binaries. Finally, the lower panels of Figure [T] show the 
relations between Mcmo and galaxy properties in the 
NSC erosion model, for nucleated galaxies with a MBH. 
Clearly, the choice of plotting the total Mcmo mass sig- 
nihcantly reduces the broadening that characterized the 
Mnsc relations (upper panels), but the scatter is slightly 
larger than in the “NSC preservation” model (middle 
panels). This is expected since Mej only equals the final 
MBH mass to first approximation (see equation [1]). 

The ancestry of a NSC includes complicated processes 
leading to its growth, disruption, and regeneration if new 
material is accreted after a disruptive merger. In Figure [5] 
we illustrate the typical evolution of NSCs as predicted 
by our fiducial “NSC erosion” model. Since the effect 
of MBH binaries and mergers on NSC’s is most impor¬ 
tant for large galaxies (where they cause the broaden¬ 
ing/bending of the correlations between NSC and host- 
galaxy properties), we choose a galaxy with spheroid 
mass ^ 2 X IO^^Mq, at which the NSC erosion becomes 
significant. In this example, the MBH grows by short¬ 
lived accretion events triggered by bar instabilities of the 
host’s galactic disk (at z ^ 6 — 7) or by major galactic 
mergers (at z ~ 1 — 3; note the corresponding growth 
of the reservoir following the merger-driven infall of cold 
gas to the nuclear regions). The NSC grows more gradu¬ 
ally through a combination of in-situ star formation (es¬ 
pecially at z ~ I — 3, when the growth of the reservoir 
triggers nuclear star formation), and infall of stellar clus¬ 
ters. The latter channel dominates over the former at 
high redshift, and is also enhanced at z ^ 1 — 3, since 
major galaxy mergers also cause violent bursts of star 
formation, which in turn enhance stellar-cluster forma¬ 
tion. Finally, at z < 1 — 2 the NSC gets eroded by a 
series of minor MBH mergers, but re-forms thanks to a 
steady infall of stellar clusters, which form as a result 
of quiescent star-formation activity as the galactic disk 
re-grows. 

Figure [3] shows the redshift evolution of the M — a 
relations for both MBHs and NSCs, in the “NSC ero¬ 
sion” model (left panels), and in the “NSC preserva¬ 
tion” model (right panels). In both models, the MBH 
growth appears inefficient at high redshifts relative to 
NSC growth. At redshift z < 2, after the epoch of bright 
quasars, MBHs become the dominant nuclear component 


in the most massive galaxies. Therefore, it is only at rela¬ 
tively late cosmic epochs that MBH mergers become effi¬ 
cient at scouring their host clusters, in the “NSC erosion” 
model. As can be seen, the “NSC preservation” model 
clearly over-predicts the mass of NSCs in the brightest 
nucleated spheroids, while the “NSC erosion” model at 
z = 0 is in good agreement with the local M — a relation 
for both MBHs and NSCs. 

4. CONCLUSIONS 

A comparison of our results to observational data re¬ 
veals that the NSC scaling correlations in the local Uni¬ 
verse carry an imprint of the merger and growth history 
of their companion MBHs. NSCs are significantly eroded 
by MBH binaries, causing a broadening of the observed 
NSC empirical correlations in high-mass galaxies and at 
low redshifts. 

The NSC scaling correlations can potentially be used 
to probe different NSC evolutionary models, and also 
to place constraints on the merger and growth history 
of their host galaxies. However, the slope of the ob¬ 
served correlations involving NSC masses is lowered by 
the inclusion, at the high-cr end, of NSCs that were partly 
eroded by MBHs. Hence, such relations cannot be used 
to put reliable constraints on different NSC formation 
models without also taking into account the scouring ef¬ 
fect of MBH binaries. 

We argue that the stellar mass removed from the center 
during mergers is of order of Mmbh, and indeed we find 
that by replacing Mnsc with Mcmo = Mnsc + Mmbh, 
the observed scaling correlations agree remarkably well 
with those predicted by models that do not account for 
the scouring effect of MBH binaries. Our findings provide 
a natural explanation to why relations between Mcmo 
and galaxy properties for nucleated galaxies show signif¬ 
icantly less scatter than relations involving Mnsc alone. 
The relations between Mcmo and galaxy properties are 
much shallower (~ a^) than the same correlations for 
MBHs (^ cr®). Our results also explain the well known 
transition from MBH- to NSC- dominated galaxies as 
one proceeds from dwarfs to giant ellipticals, without 
the need of invoking competitive feedback processes from 
young NSC’s and/or AGN activity. Therefore, the for¬ 
mation of MBHs and NSCs appear not to be mutually 
exclusive, with NSCs being modified after their forma¬ 
tion by the merger evolution of their companion MBHs. 
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